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We investigate the crossover from the Bardeen-Cooper-Schrieffer (BCS) state of weakly-bound 
Cooper pairs to the Bose-Einstein Condensate (BEC) of strongly-bound molecular dimers in a 
gas of ultracold atoms loaded on a two-dimensional optical lattice. By using the the mean-held 
BCS equations of the emerging Hubbard model and the concept of off-diagonal-long-range-order for 
fermions we calculate analytically and numerically the pair binding energy, the energy gap and the 
condensate fraction of Cooper pairs as a function of interaction strength and Hlling fractor of atoms 
in the lattice at zero temperature. 
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I. INTRODUCTION 

Several experimental groups [H-Q have observed in ul- 
tracold alkali-metal atoms the predicted 0-lll crossover 
from the Bardeen-Cooper-Schrieffer (BCS) state of 
weakly bound Fermi pairs to the Bose-Einstein conden- 
sate (BEC) of molecular dimers. In two i, il of these 
experiments the condensate fraction of Cooper pairs 
has been studied with two hyperfine component Fermi 
vapours of ^Li atoms. The experimental data of the 
condensate fraction, which is directly related to the off- 
diagonal- long-range order of the two-body density matrix 
of fermions |l2|, are in quite good agreement with 
mean-field theoretical predictions [ij, [ij] and Monte- 
Carlo simulations at zero temperature, while at finite 
temperature beyond-mean-field corrections are needed 
p^ . Recently the condensate fraction in the BCS-BEC 
crossover has been theoretically investigated for a two- 
dimensional (2D) uniform Fermi gas [17|, for a uniform 
three-spin-component Fermi gas with SU(3) symmetry 
[isj . for a 2D uniform two-component Fermi gas with 
Rashba spin-orbit coupling [19|, [20|] , and also for neutron 
matter [21|. Two years ago 2D degenerate Fermi gases 
have been experimentally realized for ultra-cold atoms in 
a highly anisotropic disk-shaped potential 
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Motivated by these recent theoretical and experimental 
achievements, in the present paper we analyze the con- 
densate fraction in the BCS-BEC crossover for a quasi- 
2D two-component Fermi gas under optical confinement, 
which gives rise to a two-dimensional square lattice [23 |. 
In particular we study the energy gap and the condensate 
fraction of Cooper pairs as a function of the interaction 
strength (or equivalently as a function of binding energy 
of pairs) and filling factor of atoms in the lattice by using 
the concept of off-diagonal-long-range-order |10l - [l3 | and 
solving the mean- field BCS equations [23|. The paper is 
organized as follows. In Section II we introduce the model 
Hamiltonian which describes two-spin-component Fermi 
atoms loaded onto a quasi-2D optical lattice. In Section 
III we discuss and solve the zero-temperature mean-field 
BCS equations as a function of the adimensional ratio be- 



tween the interaction energy per site and the tunneling 
energy, calculating the binding energy of atomic pairs, 
the chemical potential and the energy gap order param- 
eter. In particular, we compare the numerical results 
obtained by using the exact density of states with the 
analytical ones derived from an approximated density of 
states. In Section IV we calculate the condensate frac- 
tion of atomic pairs investigating the dependence of the 
condensate fraction on the relevant parameters of the sys- 
tem: scaled inter-atomic strength and filling factor. The 
paper is concluded by Section V. 



II. FERMI ATOMS ON A QUASI-2D LATTICE 

The Hamiltonian of a confined dilute and ultracold gas 
of two-component Fermi atoms is given by 
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where ipai^) is the fermionic field operator that destroys 
an atom of pseudo-spin a {a =t, J.) at the position r 
and g = Airh^as/ra is the interaction strength of the con- 
tact inter-particle potential with the s-wave scattering 
length. The external optical potential 



(2) 



produces a harmonic confinement along the z axis and a 
periodic potential 
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in the (x, y) plane, with A the wavelength of the laser light 
which determines the optical lattice [IS]- The minima 
of the lattice potential form a two-dimensional square 
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lattice with sites in the positions Rj = ai = a{ix,iy), 
where a — A/2 is the lattice spacing and i — {ix,iy) is a 
2D vector of integer numbers. 

Using the set of Wannier functions in the lowest Bloch 
band |24], where the Wannier function Wi{x,y) is max- 
imally localized at site Ri, we can expand the fermionic 
field operator as: 
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where cia and cia obey the usual Fermi anti-commutation 
relations, and = y^h/{mujz) is the characteristic 
length of the strong harmonic confinement along the z 
axis, which induces a quasi-2D confinement if huJz is much 
larger than the other energies of the system. Under these 
conditions, the Hamiltonian ([T]) can be written as 
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where (ij) means nearest neighbor sites. 



t = -J dxdyW-:{x,y) 
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is the hopping parameter {t > 0), i.e. the tunneling en- 
ergy between nearest neighbor sites, and 



U=-^ f dxdym{x,y)\' 
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is the on-site strength of the inter-atomic interaction. 
fiio- = cjl^Cio- is the number operator which describes the 
number of atoms with spin a at the site i, and conse- 
quently the total number operator reads 
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Notice that Eq. (7) holds under the conditions \as \ <^ Uz 
and \as\ <C a, which ensure the absence of confinement 
induced resonance (23| and no distorsion of Cooper pairs 
due to neighbor valleys of the optical confinement. In the 
Hubbard-like Hamiltonian ([5]) we have not included the 
tunneling energies between sites which are not nearest 
neighbor because they are exponentially suppressed. We 
have also assumed the on-site one-body energies to be the 
same on all sites and therfore dropped them as irrelevant 

M. 



III. MEAN-FIELD BCS EQUATIONS 

It is well known that the BCS state appears only in 
the case of an attractive strength, i.e. U < In 
the past the negative-C/ Hubbard Hamiltonian has been 
investigated by various authors [1^ as a model for high- 
Tc superconductivity. More recently, it has been used 



to study the BCS-BEC crossover on 2D and 3D lattices 
both at zero [26,n28i] and finite temperature [1^ [s^l . As 
stressed in the introduction, motivated by recent theoret- 
ical and experimental achievement with ultracold atoms 
in optical lattices, here we reconsider the 2D negative-C/ 
Hubbard Hamiltonian to investigate the pair condensa- 
tion, and in particular the condensate fraction of Fermi 
atoms in the 2D lattice at zero temperature. Note that 
he condensate fraction has been calculated by Kujawa 
(27} in the 3D square lattice with a generalized Hubbard 
model, but only in the special case \U\/t = 00 . In the fol- 
lowing sections we calculate, as a function of |C/|/t and of 
the filling factor, the energy gap and condensate fraction 
in the 2D square lattice, analyzing also the pair bind- 
ing energy, which is always finite in the 2D BCS-BEC 
crossover. 

We start by decoupling the interaction Hamiltonian of 
Eq. ([5]) in both normal and anomalous channels 31] 
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We also assume 
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and introduce the (real) mean- field, site-independent, 
gap order parameter 

A = -U{c+c+) = -(7(ci^cit) . (11) 

In this way we obtain the mean-field Hamiltonian 
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where Ng is the number of lattice sites. 

In the dual space of wavevectors k = {kx, ky), setting 
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where Cko- destroys an atom of spin <t and wavevector k, 
the mean-field Hamiltonian (|13p becomes 
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where 



ek = — 2i (cos (fc^a) + cos [kyo)) 
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is the single-particle energy. We stress that we are con- 
sidering only the lowest Bloch band. This single-band ap- 
proximation for the BCS-BEC crossover is reliable since 
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the crossover occours at raagnetic fields that are rela- 
tively far away from the Feshbach resonance underlying 
it [32| . Moreover the app roximation is reliable under the 
following conditions [32], [s^ : i) there are no more than 
two fermions per site; ii) the two lowest bands do not 
overlap, implying that Vq 3> Er, which means 8i <C Er, 
and \U\ <^ Eg. Er = h^k\/ (2to) is the recoil energy with 
Ul = In I a the wavevector of the 2D optical lattice, and 
Eg is the energy gap between the first and the second 
Bloch band. 

We calculate the thermodynamic potential 
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where fi is the chemical potential which determines the 
average number N — (N) of fermions, by introducing the 
Bogoliubov canonical transformation: 

"k = Uk Ckt " c^k^ , /3k = ?^kC-ki + WkCk| , (17) 



= 1- After the 



where Wk and Uk sue real and 
mimimization of f2 with respect to /x and A we recover the 
standard BCS equation [2J, [3l| for the average number 
of particles per site 



= V^k , 

k 

and the familiar BCS gap equation 
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where the quasi-particle amplitudes Uk and Wk are given 
by 



1 



1 



Ck - h 
Sk 



(20) 



and 



where 
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^. Here the Bogoliubov energy reads 
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is the effective chemical potential which takes into ac- 
count the Hartree interaction. The effective chemical po- 
tential h and the gap energy A are obtained by solv- 
ing equations (1181) and (jl9p . In the continuum limit 

Ek- 

of states (DOS) per site 
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FIG. 1: (Color online). Scaled binding energy Esji as a 
function of the scaled interaction strength 1(71/4, with t the 
tunneling energy. Solid lines are the results obtained with the 
exact density of states (exact DOS) given by Eq. while 
dashed lines are the results obtained with the approximate 
density of states (approx DOS) given by Eq. H27|) . 



where BZ — [— 7r/2,7r/a] x [— 7r/a, ■7r/2] is the first Bril- 
louin zone, K(x) is the complete Elliptic integral of the 
first kind and Q{x) is the step function, the number equa- 
tion ((T8)) and the gap equation ((T9| can be written as 
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As discussed in ^25'] , quite generally in two dimensions 
a bound-state energy Eb exists for any value of the neg- 
ative interaction strength I] . For the contact potential 
the bound-state equation in the lattice is 
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where eo — —At is the lower value of the single-particle 
energy ek, occurring at k = 0. If we approximate the 
true DOS with a constant value in the interval [— 4t,4t], 
i.e. 
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that ensures the normalization 
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the bound-state equation can be solved analytically giv- 



ing 
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FIG. 2: (Color online). Scaled eflFective chemical potential ft/t 
(upper panel) and scaled chemical potential n/t (lower panel) 
as a function of the scaled interaction strength \U\/t, with t 
the tunneling energy. Results obtained for three values of the 
filling factor x = n/2. Filled circles are the results obtained 
with the exact density of states given by Eq. (|23p . while 
lines are the results obtained with the approximate density of 
states given by Eq. (f27)l . 



In Fig. [T]we plot the binding energy EB/t as a func- 
tion of the interaction strength \U\ obtained with this ap- 
proximate formula (dashed line) . For comparison we plot 
also the exact result (solid line), obtained by numerically 
solving Eq. (1^51) . The figure shows that the agreement 
between the two curves is extremely good. The BCS- 
BEC crossover is governed by the adimensional param- 
eter \U\/t or equivalently by the scaled binding energy 
Eb /t. The limit of large tunneling and small interaction 
\U\/t ^ 1 corresponds to the ECS regime where EB/t 
is close to zero. Instead the limit of strong localization 
and large interaction \U\/t 1 corresponds to the BEC 
regime where EB/t is large. 

The quite good agreement between the solid curve and 
the dashed curve of Fig. [l]suggests that one could use the 
approximate DOS to study various ground-state proper- 
ties of the system in the BCS-BEC. Within the approxi- 
mation of a constant DOS in the band, i.e. Eq. ((27)) . the 
number density equation and the gap equation read 
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It is then straightforward to plot (see Fig. [2]) the effec- 
tive chemical potential h (upper panel) and the chemical 
potential fi (lower panel) as a function of the scaled in- 
teraction strength \U\/t, for different values of the filling 
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FIG. 3: (Color online). Upper panel: Scaled energy gap A/t 
as a function of scaled interaction strength \U\/t with t the 
tunneling energy. The three curves correspond to five different 
values of the filling factor x = n/2. Lower panel: Scaled 
energy gap A/t as a function of filling factor x — n/2, where 
the three curves correspond to three different values of the 
scaled interaction strength \U\/t, with t the tunneling energy. 
Filled circles are the results obtained with the exact density of 
states given by Eq. (|23|l . while lines are the results obtained 
with the approximate density of states given by Eq. (p7)) . 



factor X — n/2{0<x<l). In the figure the lines are 
obtained by using Eqs. ([50]) and ([50)1 based on the ap- 
proximate DOS of Eq. (P7)) . while the filled circles are 
obtained by using Eqs. (p4|) and p5|) with the exact DOS 
ofEq. (1121). 

Fig. [2] shows that at half filling (x = 0.5) the effec- 
tive chemical potential h remains always constant and 
equal to zero, and the corresponding chemical potential 
/i follows the simple law /i = —\U\/2. Moreover, the 
lower panel of Fig. [2] shows that, at fixed filling factor x, 
the chemical potential /i as a function of U is close to a 
straight line (it is true straight line only for x — 0.5) and 
approaches ^ ~ —\U\/2 for large \U\. 

In Fig. [3] we plot the energy gap A vs interaction 
strength \U\ (upper panel) and vs filling factor x (lower 
panel). The upper panel shows that, at fixed filling fac- 
tor X, the energy gap A grows by increasing the scaled 
interaction strength \U\/t, that is by increasing the lo- 
calization. Instead, the lower panel shows that, at fixed 
scaled interaction strength |C/|/t, the scaled energy gap 
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A/t reaches its maximum value at half filling x = 1/2, 
i.e. when on the average there is one atom per site. This 
effect is clearly seen in the lower panel of Fig. [3] where we 
consider three values of Notice that the behavior of 

A as a function of x is perfectly symmetric with respect 
to X = 1/2 (half filling). Also in Fig. [3] the agreement 
between the results obtained with the exact DOS and the 
ones calculated with the approximate DOS is quite good, 
and it improves by increasing \U\/t. Motivated by this 
finding, in the remaining part of the paper we use the 
approximate DOS, which is much simpler for numerical 
computations and produces analytical results. 

IV. CONDENSATE FRACTION 

The main task of the paper is to analyze the conden- 
sate fraction of fermions. As shown by Yang \Td], the 
BCS state guarantees the off-diagonal-long-range-order 
prj of the Fermi gas, namely that, in the limit wherein 
both unprimed coordinates approach an infinite distance 
from the primed coordinates, the two-body density ma- 
trix factorizes as follows: 

(^+(rl)v3+(r'2)V';(ri)^t(r2)) (32) 
= (^+(r;)^+(r'2))(^,(ri)V't(r2)). 

The largest eigenvalue of the two-body density matrix 
(|32]) gives the number of pairs in the lowest two-particle 
state, i.e. the condensate number of Fermi pairs 
[T2j |. In this way, the number No of condensed fermions 
is given by 

ij 

(33) 

Notice that, as said above, A'o counts the number of con- 
densed fermions: < Nq < N [31 and not of condensed 
pairs. It is then quite easy to show that the the conden- 
sate number of atoms per site is 

^o = 2l^E"k«k- (34) 

k 

With the help of Eqs. (|20)) this number is thus given by 
A2 Z-^* 1 

and using the approximate DOS of Eq. (f27|) it reads: 

A / f4t~h\ fAt + h\\ 

^ IW 1^''*^" ) + ) ) ■ 

(36) 

Fig. 21 shows the condensate fraction no/n of fermions, 
calculated with Eqs. ([50)) . (pij) and ([55]) . as a function of 
scaled interaction strength \U\/t for three values of the 
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FIG. 4; (Color online). Condensate fraction no/n as a func- 
tion of scaled interaction strength \U\/t with t the tunneling 
energy. The three curves correspond to three different values 
of the filling factor x = n/2. 



filling factor x. We have verified that the plotted results 
are in good agreement with the ones obtained by using 
the exact DOS, except in the case of very small values 
of |C/|/i. In any case, the condensate fraction no/n van- 
ishes when the scaled interaction strength \U\/t goes to 
zero. Moreover, as shown in the figure, the condensed 
fraction grows very fast for values of the scaled interac- 
tion strength \U\/t < 8, it shows a shoulder, and then it 
reaches its asymptotic value ng/n ~ 1 — x rather slowly. 

This result is confirmed in the upper panel of Fig. 
[5l where we report the condensate fraction no/n as a 
function of the filling factor x at fixed scaled interac- 
tion strength \U\/t. The figure clearly shows that no/n 
ranges from one to zero, being extremely close to one for 
X <^ 1 and approaching zero as x goes to 1. This means 
that there is a full BEC-BCS crossover by increasing x 
at constant scaled interaction strength |?7|/t. Moreover, 
if the scaled interaction strength \U\/t is large, the con- 
densate fraction no/n follows a straight line during the 
BEC-BCS crossover. For the sake of completeness, in 
the lower panel of Fig. [5] we plot also the number no of 
condensed atoms per site as a function of the filling fac- 
tor X. The results show that the curves of no vs x have 
a behavior similar to those of A vs a; (see Fig. [3]). In 
the limit \U\/t — > 00 one finds that no = {1 — x)2x and 
consequently no/n = {1 — x). 

Finally, we observe that, after a simple rescaling of the 
chemical potential, namely h ~ h+4t, in the limit t — ?► 00 
with ta^ — > TTh?/m, Eq. becomes the condensate 

number equation found in Ref. [17i] for the 2D uniform 
superfluid Fermi gas. 



V. CONCLUSIONS 

By using the mean-field extended BCS theory and the 
concept of off-diagonal long-range order, that is the ex- 
istence of a macroscopic eigenvalue of the two-body den- 
sity matrix, we have investigated the condensate fraction 
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FIG. 5: (Color online). Upper panel: Condensate fraction 
no/n as a function of the filling factor x — n/2. Lower 
panel: Number no of condensed atoms per site. The curves 
correspond to three different values of the scaled interaction 
strength with t the tunneling energy. 



of fermionic pairs in a uniform 2D Fermi gas. We have 
shown that the condensate number uq of fermi atoms 
per site is extremely useful to characterize the BCS-BEC 



crossover, that is induced by changing the adimensional 
ratio \U\/t between the interaction energy \U\ and the 
tunneling energy t. In particular, we have found that 
both the scaled binding energy Eb /t of atomic pairs and 
the condensate fraction no /n grow by increasing the ratio 
\U\/t at fixed filling factor x — n/2 (with n the average 
number of fermions per site). In addition, our results sug- 
gest that fixing the ratio \U\/t, or equivalently the scaled 
binding energy Es/t, there is a full BEC-BCS crossover 
by increasing the filling factor from zero to one. Finally, 
we have found that the analytical results obtained by us- 
ing an approximate density of states are in quite good 
agreement with the numerical ones deduced from the ex- 
act density of states. In our calculations we have used 
the mean-field theory and it is important to stress that 
recent Monte Carlo simulations have shown that, at zero- 
temperature, beyond-mean-field effects are negligible in 
the BCS side of the BCS-BEC crossover while they be- 
come relevant in the deep BEC side 0, [s^l • In any case, 
we think that our mean-field results, and the reliable an- 
alytical formulas we have obtained, can be of interest for 
near future experiments with degenerate gases made of 
alkali-metal atoms confined in quasi-2D optical lattices. 

In this paper we have investigated zero-temperature 
pair condensation. According to the Mermin- Wagner 
theorem [2^, for an infinite 2D system there is conden- 
sation (ofF-diagonal-long-range order) only at zero tem- 
perature. However, for a finite 2D system condensation 
could be possible also at non-zero temperature. The in- 
vestigation of this issue, which requires a beyond mean- 
field approach, for 2D fermions in a lattice is in progress. 
Another puzzling issue is the filling of the second Bloch 
band: we plan to investigate its effect on pair conden- 
sation by analyzing a multi-band version of the present 
theory. 
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